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VQ , We investigate theoretically the fluid-crystal coexistence of solutions of globular charged nanopar- 

f^ ' tides like proteins and inorganic colloids. The thermodynamic properties of the fluid phase are com- 

f^ ' puted via the optimized Baxter model. This is done specifically for lysozyme and silicotungstates 

^S| . for which the bare adhesion parameters are evaluated via the experimental second virial coefficients. 

The electrostatic free energy of the crystal is approximated by supposing the cavities in the inter- 
O .' stitial phase between the particles are spherical in form. In the salt-free case a Poisson-Boltzmann 

equation is solved to calculate the effective charge on a particle and a Donnan approximation is 
used to derive the chemical potential and osmotic pressure in the presence of salt. The coexistence 
data of lysozyme and silicotungstates are analyzed within this scheme, especially with regard to 
the ionic-strength dependence of the chemical potentials. The latter agree within the two phases 
provided some upward adjustment of the effective charge is allowed for. 
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I. INTRODUCTION 



' O ■ One current view of protein crystallization centers on the second virial coefficient B2 being a relevant quantity 

^ ' determining the onset of crystallization [H Q, y, |j] . There exists a crystallization slot of negative B2 values vifhich 
expresses a necessary range of solution conditions for adequate crystals to grow. A negative value of B2 implies a 
Baxter stickiness parameter as it is conventionally defined and we here denote by tq. Thus, in a similar vein, tq has 
been correlated with the solubility of nanoparticles in explaining fluid-crystal coexistence curves |^ Q . 

The free energy of a suspension of particles cannot, of course, depend on B2 alone. In a recent paper [7| we 
introduced a new analytical theory for protein solutions in which the real fluid is replaced by a suspension of spheres 
with an appropriately chosen adhesion of the Baxter type. The stickiness parameter r in the latter is computed by 
^Sl \ a variational principle for the free energy. In our optimized Baxter model, r is not at all identical to tq but depends 
^^ not only on the ionic strength but also on the protein concentration. In Ref. |3, Rosenbaum et al. plotted tq 

^^ logarithmically as a function of the nanoparticle concentration which effectively coarse-grains the experimental data 
^C , they show. If we zoom in on their curve, there is a lot of fine detail which we here argue to be related to the fact 
that T is a better similarity parameter. In particular, we seek to understand the ionic-strength dependence of the 
fluid-crystal coexistence curves by going beyond theory based solely on tq. 

We have recently tested the optimized Baxter model on a system of spheres interacting via an attractive Yukawa 
potential analyzed by computer simulations '8|. The stickiness parameter r, evaluated by optimizing the free energy, is 
i-pj ■ indeed a useful similarity variable for gaining insight into the pressures and chemical potentials from the simulations. 
(~| ' The magnitudes of these quantities are also well predicted by the optimized Baxter model. 

O ■ However, here we will not focus on the variable t and the fluid phase but rather on the coexistence itself. The 

^ ' systems we study are assumed to have a short enough range so that the coexistence between two fluid phases is 
circumvented. An a priori theory is problematic because we would need a quantitative theory of the crystal phase in 
terms of postulated attractive forces which are currently unknown. Theoretical efforts exist in the literature ^, 10. ll| 
at the expense of introducing unknown parameters which we want to avoid here, 
jrt ' In practice, it may be very difficult to achieve ideal thermodynamic equilibrium between the liquid phase and some 

■ ■ ' crystalline state. Equilibrium may not have been reached, the crystal could be heterogeneous and the formation of 
aggregates could complicate the attainment of equilibrium (see, for instance, the discussion by Cacioppo and Pusey 
on lysozyme |l_25). Nevertheless, it may still be useful to assume equilibrium is ideally attained provided our goal 
is sufficiently modest. The balance of chemical potentials has been used before to acquire information about the 
crystal from the solubility in the fluid phase |3- Our concern here will be to try to gain insight into the ionic- 
strength dependence of the thermodynamic properties of the crystal. We may argue that this dependence could 
be approximated by a Donnan equilibrium so it would not be very sensitive to the precise crystal habit adopted. 
We therefore compute the protein chemical potential and osmotic pressure of the coexisting liquid phase at the 
experimentally determined solubility with the help of the optimized Baxter model. We then investigate whether their 
dependence on the electrolyte concentration agrees with that predicted by a simple crystal model. 



II. OPTIMIZED BAXTER MODEL 

We first discuss liow we obtain the bare adhesion parameters via the second virial coefficient, and then summarize 
the optimized Baxter model which is a liquid state theory at finite concentrations Q ■ We consider a system of charged 
nanometer-sized particles (e.g. proteins or nanocoUoids) in water with added monovalent salt of ionic strength /. We 
suppose the particles are spherical with radius a. The charge is distributed uniformly on the particle's surface. 
For convenience, all distances in this section will be scaled by the radius a and all energies by ksT, where fc^ is 
Boltzmann's constant and T is the temperature. Because monovalent ions (counterions and salt ions) are present in 
solution, the Coulomb repulsion between the particles will be screened and it is here given by a far-field Debye-Hiickel 
potential [Tl . The effective number Z^ff of charges on the sphere (taken to be positive) will here be computed in the 
Poisson-Boltzmann approximation. We let the attraction between two particles be of range much shorter than their 
radius, and we model it by a potential well of depth Ua and width (5 <C 1. The total interaction U{x) between two 
particles whose centers of mass are separated by an actual distance r is thus of the form 

(X) <x <2 

U{x)^{Udh{x)-Ua 2<x<2 + 5 , (1) 

Udh{x) x>2 + 5 

r 

x= -, 
a 



with Debye-Hiickel interaction 



g-w(x-2) 



Udh{x) = 2C . (2) 
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Here, ^ = ^ I -rr^ ) and oj = Ka, which are given in terms of the Debye length k ^ defined by k^ — SttQI and 

the Bjerrum length Q = q^/eksT, which equals 0.71 nm in water at 298 K (e is the permittivity of water, q is the 
elementary charge); w = 3.28av/, if the radius a is given in nm and the ionic strength / in M. We suppose 1-1 
electrolyte has been added in excess so / is the concentration of added salt. We have derived the effective charge 
qZeff in the Poisson-Boltzmann approximation [2| 

Here Ei{x) is the exponential integral defined by -^1(2:) = / dtt^^e^*^ and qZ is the actual charge per particle. 

We suppose that the bare charge on the particles as a function of the ionic strength is known from experiment, so 
the only unknown parameters are Ua and 5 which are chosen to be independent of /. The latter are determined by 
fitting preferably complete experimental data of the second virial coefficient B2 as a function of the ionic strength / 
at constant pH to B2 computed numerically with the help of the expression 



B2 = 2TTa / X dx 



(1 - e-^(-)) (4) 



using Eq. QJ. We have previously done this for a wide variety of B2 data on lysozyme at two values of the pH (4.5 
and 7.5) and we were able to obtain very good fits Q (see e.g. Fig. 1 which is discussed in section HI.B). 

It is important to stress that though there are two adjustable parameters 5 and Ua, the actual fit in practice depends 
almost solely on adjusting the single combination SexpUA- This is because a convenient analytical approximation of 
the second virial turns out to have the form [3 

^2 .,i + ^_V«<5e^- (5) 



B"S - ' 2lo 



'2 



and is able to describe the experimental data on lysozyme quite well with an appropriate value of JexpC/yi. Here, 
B^^ is the second virial coefficient pertaining to hard spheres. We note that Eq. disagrees starkly with an 
approximation put forward earlier [13, both with regard to the pure electrostatic and the adhesive contributions. In 
particular, the third i.e. adhesion term in Eq. ((S)) is not at all independent of the ionic strength but rather diminishes 
fast as the electrolyte concentration is lowered. Furthermore, the pure electrostatic term cannot be derived from a 
Donnan equilibrium as we point out in Section IV. 



At high salt concentrations, the parameter ^ becomes small owing to screening so B2 becomes lower than the hard 
sphere value as can be seen from Eq. ((SJ. Nevertheless, the electrostatic repulsion still exerts itself, so an effective 
adhesion parameter we may wish to introduce would be smaller than the bare value. We therefore adopt a similar 
strategy to the liquid state at finite concentrations by first introducing a suitable reference state amenable to analytical 
computation Jj . This is a solution of hard spheres whose radius is still a but with a Baxter adhesion potential whose 
strength is defined by a suitable stickiness parameter t. The statistical properties of this suspension as a function of 
the volume fraction of spheres rj (= Ana^ /3 times number density) may be solved in the Percus-Yevick approximation 
|l4| . The parameter r is adjustable and is computed via a variational principle for the free energy. The latter may 
be written as a functional expansion in terms of the so-called blip function which is the difference in Mayer functions 
of the respective interactions (Eq. ^ and the Baxter interaction) |3, U^ ■ We set the first-order deviation from the 
free energy pertaining to the reference state equal to zero. This determines r which depends not only on the well 
parameters S and Ua and electrostatic variables lu and ^ but also on the volume fraction of nanospheres. It is given 

by a 



T 



.Ua^~t+tt2^ -Ij (l + (l + H)e) + (e^^e~«-l) , (6) 

where 

e^6-K-^ [(1 + dH)Pi + HP2] , (7) 

8 , ,, 16 / M \ 

8 16 

P2 = ^(2 + c^<5)M+^ln(l + M), (9) 

M = Ce-"V4, (10) 

Ki = 2 (e^'^e'^'^'''^' - l) {1 + {1 + H) S) + 2 (e^^e'^^''"^^ - iVl + (1 + H) e) , (11) 

2t(1-77) V 12 12 l-V 2(l-r;)2Ay ^ ' 



and A is given by 



v/2 1 V , V_^^ (,3) 



(1-77)2 A 1-T] 12 



Note that r is readily obtained by iteration. One starts with initial values for t and e and then calculates A, H and 
Ki from Eqs. (|11|I - H13|I . Then, a new value of e at fixed H is computed iteratively with the help of Eqs. and 
(|ll|l . Next, a new value of r is given by Eq. © and then the cycle is repeated until the variables become stationary. 
Having obtained the effective adhesion parameter t, we simply calculate thermodynamic properties of the reference 
state within the Percus-Yevick approximation. The free energy of the actual system does deviate slightly from that 
of the reference state but we have shown that the deviations are very small Q . To compute the osmotic pressure 11 
we use the result from the compressibility route [lj| 

n^;o ^ Tj(l + n + Tj^) 7^2(1 -fry/2) r;3 

keT (1 - 77)3 (1 - 77)2 '' ^ 36 ■ ^ ' 

When the roots of Eq. (jl3|l are complex, the pressure cannot be determined for the physical realization of the liquid 
state breaks down, at least with the Percus-Yevick approximation. The chemical potential /i of the spherical particles 
is determined by using the pressure from Eq. (|14ll and the Gibbs-Duhem equation at constant temperature [l6j 



Here 



J= 3 . _ 3,(1 + 4,) 6,(2 + ,) _ ^8^ _ 6(r - r.)^ 

2' (1-,) (l-r/)2 1-, r,(l-6T,) 

, 6Te(18rT, - 1)2 
1 - 6t. 



A(l - 77) - T-i 



A(l - 77) - 18t, 



r-i - 18r, 



(16) 



is the contribution to the chemical potential that vanishes in the hard-sphere limit (r — > 00) and 



where h is Planck's constant and m is the mass of a sphere. The critical value of r below which there is a range of 
densities where there is no real solution of A, is given by 

r. ^ ^^ (18) 



III. SOLUBILITY CURVES: CHEMICAL POTENTIAL OF THE FLUID PHASES 

A. Method 

Since we suppose the crystal is in thermodynamic equilibrium with the fluid, the protein chemical potentials as well 
as the osmotic pressures in both phases are uniform. The chemical potential of the counter and co-ions must also be 
uniform but we will address this issue later within a Donnan equilibrium. Solubility data from experiment represent 
the particle concentration in the fluid phase as a function of the pH and the salt concentration. Thus we compute 
the chemical potential and the osmotic pressure of the solution with the help of the optimized Baxter model of the 
previous section. We have done this in two cases of nanoparticles where we have sufficient experimental data on the 
second virial coefficient to evaluate the well parameters Ua and 5 with sufficient accuracy. 

B. Lysozyme 

The protein hen-egg-white lysozyme has been well characterized in aqueous solutions of simple electrolytes. We 
here choose the effective radius a such that the volume of the model sphere is equal to the volume of a lysozyme 
molecule in the tetragonal crystal state. The latter is determined from the water content of the tetragonal crystal 
(0.335 mass fraction T7]), the crystal volume per protein molecule (29.7 nm'^ [l3)i the density of the crystal (1.242 
10'^ kg m"'' [T3|) and the density of water (0.998 10"^ kg m^'^). Thus we have a = 1.61 nm and note that this is about 
0.1 nm less than the value of 1.7 nm we used previously 0, which was based on approximating the protein by an 
ellipsoid of dimensions 4.5 x 3.0 x 3.0 nm [13. For the sake of consistency we here use the single value a = 1.61 nm 
in computations pertaining to both phases. 

The experimental data for the second virial coefficient of lysozyme have been discussed by us at length previously 
Q and are presented in Fig. 1. For details on determining the parameters Ua and S of the attractive potential we 
also refer to Ref. Q. Since we are using a smaller effective radius here, we deduce the values Ua = 2.89 and 5 = 0.182 
which are somewhat different from those derived earlier 7] . The values of the bare charge qZ of a lysozyme molecule 
as a function of the ionic strength are the same as those used in Ref. l7| i.e. they are determined by interpolation 
from hydrogen- ion titration data in KCl 2Q|. We assume that KCl and NaCl (see below) behave identically in an 
electrostatic sense. The effective charge does differ slightly because it is a function of a (see Eq. Q). We again use 
the lowered effective charge Z — Z^ff — 1 instead of the effective charge Z^ff in order to fit B2 accurately at lower 
ionic strengths when it is dominated by electrostatics. We set Ua and 6 to be independent of the pH. 

Accurate data on the solubility S" as a function of the NaCl concentration have been obtained by Cacioppo and 
Pusey [13 using column beds of tetragonal microcrystallites of lysozyme in a range of pH and temperatures. We here 
employ their data at 298 K and at three representative values of the pH. (See Tables I, II and III.) The ionic strength 
/ in M is determined from the ionic strength in %w/v by the relation /(M)= 0.06 -I- 0.171 /(%w/v). Here, the value 
0.06 accounts for the effective ionic strength of the 0.1 M sodium acetate buffer used and 0.171 = 10/MNaCi where 
Mj^aCi = 58.44 g mol~^ is the molar mass of NaCl. The dimensionless parameter uj is then given hy uj — 5.28\/7, 
where / is given in M, and ^ = 0.220(Z/(1 + u;))^. 



The volume fraction 77 of protein in the hquid phase is given by 77 = SN avvi/M , where Na is Avogadro's number, 
Vq = 47ra'^/3 is the volume of a lysozyme molecule and M = 14.3 kg mol~^ [2l| is the molar mass of lysozyme. The 
parameter r describing the effective adhesion is determined as described in section II (see Eqs. (|H|>- (|1^ 1 using the 
values Ua — 2.89 and 5 — 0.182. Then, the dimensionless chemical potential {fi — ijLo)/kBT and the dimensionless 
osmotic pressure Hvo/ksT are determined from Eqs. 1151 and [TH respectively. (See Tables I, II and III). Fig. 2 shows 
the chemical potential as a function of the ionic strength / at three different values of the pH. Fig. 3 shows the 
osmotic pressure under the same conditions. 

C. STA 

The next system we consider is silicotungstate (STA) in water with three different kinds of added salt: NaCl, 
HCl and LiCl. STA molecules are spherical, more or less, (see Fig. 2 in Ref. |23|) with an effective diameter of 
1.1 nm |23,|23|, so we set a = 0.55 nm. The structural formula for the polyanion SiWi204(^ implies a molar mass 
MsTA = 2874.3 g mol~^. We assume that the pH is low enough for the molecule to be fully dissociated, i.e. Z = A. 

We determine the well parameters Ua and S for the attractive interaction by fitting experimental data of the second 
virial coefficient in the same way as was done for lysozyme Q, except we now do not adjust Z^ff- The second virial 
coefficients for Li4STA, H4STA and Na4STA are taken from Ref. |25j and plotted in Fig. 4. In each case, the added 
salt is XCl, where X represents the counterion of the crystal. The values of the dimensionless parameters uj = 1.80\/T, 
^ = 0.645(Ze///(l + uj))'^, Z and Zeff pertaining to the data in Fig. 4 are given in Table IV. We have set S = 0.05. 
A least-squares fit to the data represented in Fig. 4 then gives Ua = 3.30. In fact, there is a range of combinations 
of 5 and Ua that yield almost the identical curve as long as i5exp Ua ~ 1-36 and J ^ 1, so our choice oi 5 — 0.05 is a 
bit arbitrary. This similarity with respect to the sole parameter 5 exp Ua is in accord with our approximation for B2 
give by Eq. © . 

The solubilities for Li4STA, II4STA and Na4STA have been measured by Zukoski et al. [23, where the same 
electrolytes are used as in the measurements of 52- (See Tables V, VI and VII.) The volume fraction r] of STA 
is given by 77 = SNavo/Mx^sta, where S is the solubility of STA (note that here it is given in g/ml, whereas for 
lysozyme it was given in g/1), vq = 47ra^/3 is the volume of an STA molecule and Mx^sta is the molar mass, where 
X again represents the counterion in the respective cases. We have Mh^sta = 2878.3 g mol~^, M^i^sTA = 2902.0 g 
mol~^ and MMa^STA = 2966.2 g mol~^. The stickiness parameter r is determined by the method described in section 
II (see Eqs. (|S|)- (tT^ '). using the values Ua = 3.30 and 5 = 0.05. The chemical potential and the osmotic pressure are 
again determined from Eqs. 1151 and [HI respectivelv. (See Tables V, VI and VII). We display these thermodynamic 
variables as a function of the ionic strength in Figs. 5 and 6. 

IV. CRYSTAL MODEL: DONNAN EFFECT 

Having computed the thermodynamic properties of the fluid phases of lysozyme and STA, and hence those of the 
respective crystal phases under the assumption of equilibrium of the two phases, we now attempt to gain insight into 
them by introducing a simple model for the crystal. In the latter the spherical particles either touch or are very 
close. There are thus minute " surfaces of interaction" where the forces between two nearby spheres are predominantly 
attractive. It is therefore reasonable to write the thermodynamic potential J7 of a crystal of N spheres in a volume V 
as a superposition of attractive and electrostatic contributions to a first approximation 

^ = \k{c)^^^—^ + Nf,l{c,h)+n{SJ)V-^i{SJ)N (19) 

The crystal is immersed in a large reservoir at a constant osmotic pressure 11 and chemical potential /i containing a 
saturated solution of nanospheres at a solubility S and ionic strength / (11 and ji are given by Eqs. ifHI) and l(T^ 
respectively). The crystal has elastic properties denoted by the modulus k which depends on the density c — N/V 
and the crystal would have a volume Vq in the absence of electrostatic forces (|V — Vb| <C Vq). Actually, the form of 
the elastic energy is more complicated and depends on the precise crystal habit 26j but the simple harmonic form 
in Eq. H19|) suffices for our purposes. There is a Donnan equilibrium (see below) which leads to a salt concentration 
Ic within the interstitial region in the crystal. We adopt a continuum approximation: the electrostatic free energy 
Nfei is computed for a lattice of charged spheres embedded in a solvent of uniform permittivity e and electrolyte 
concentration L. 



At equilibrium, fl must be minimized {dfl/dV — 0; dil/dN — 0) so that 



, ,'V-Vo\ 1 dfc fV-VoY , , 



Idfc fV-Vo\ 



2 dc V V'c 

We have introduced the electrostatic counterparts of the osmotic pressure and the chemical potential of a charged 
sphere in the crystal phase on the right hand sides of Eqs. (|20(l and pifl . In Eq. H2Q|I the elastic term proportional 
to k may easily be of order Il^i but the quadratic form is negligible. In view of the fact that Ilg; = 0(c/ie/), we then 
have /i ~ /ie/ to a good approximation from Eq. (|21|l . In effect, as we change the ionic strength of the fluid phase, 
the solubility S and the salt concentration /^ within the crystal readjust themselves whereas the volume V remains 
virtually constant. The chemical potential is modified only by virtue of the change in electrostatic shielding about 
a sphere in the lattice. But a substantial hydrostatic pressure may be exerted within the crystal as we decrease its 
volume a bit. 

Next, we compute the electrostatic properties of the crystal. The coUigative properties of salt-free polyelectrolytes 
are often addressed in terms of a cell model in which a test cylinder is surrounded by a boundary of similar symmetry 
on which the electric field vanishes [23 . The boundary effectively replaces the effect of the surrounding particles on 
the test particle. This picture is reasonable at low volume fractions but must break down at high concentrations when 
the electric field is highly heterogeneous. In the latter case, one of us has advocated focusing on the voidlike regions 
instead of on a test particle (see Ref. p8j which deals with a hexagonal lattice of DNA at very high concentrations). 
Thus, in a crystal of spheres we may distinguish very small regions between particles that almost touch which we 
view as thin boundary layers, and larger voids which we will simply approximate by spheres. (We are here concerned 
with spheres of high charge density which leads to counterions being "condensed". At low charge densities, it is 
possible to give a more general analysis; see the Appendix). Discrete charge effects should prevail when evaluating 
the electrostatics of the boundary layers. These energies are here assumed to be independent of the ionic strength 
since the relevant scales in the boundary layers are very small in crystals of nanoparticles. 

We therefore first solve the Poisson-Boltzmann equation for a charged void or spherical cavity of radius b without 
salt and then discuss the effect of monovalent salt via a Donnan equilibrium. The charge density on the surface of 
the cavity is uniform and the total number of charges is Z. In view of electroneutrality there are Z counterions in the 
cavity, each bearing charge —q. Within a mean-field analysis, the counterion density p{r) inside the cavity is given 
by a Boltzmann distribution in terms of the electrostatic potential '^(r) at a distance r from its center 

p(r)=pe«*/'=-^. (22) 

We choose ^ = at the center of the cavity so that p is the actual charge density there. The charge density —qp is 
also related to ^ by Poisson's equation 

A* = 1^ (23) 



leading to the Poisson-Boltzmann equation |27| which we conveniently express in the scaled form 

ib"ix) + -ij'{x) = e^. (24) 



X 



Here, we have defined A ^ = AnQp, x = r/X and ip = q'^ /ksT where A may be interpreted as a screening length. 
The two additional boundary conditions are 



owing to symmetry, and 



V''(0) = (25) 



i,(i).A.f 



signifying the relation between the electric field and the charge density at the surface of the cavity. 

For small x, Eq. (|24|l admits a series expansion i!{x) = Ax^ + Bx'^ + ... with A = 1/6 and B = 0(1) independent of 
the value of the dimensionless variable A. As A tends to zero, Eq. (|25|l reduces to the condition of electroneutrality. 



Electrostatic screening vanishes in this hmit and there are no counterions "condensed" on the surface of the sphere. 
It is straightforward to solve Eq. H24I) numerically starting with ^|J{x) -^ ^x^ as x ^ 0. We have fitted the solution 
to the convenient approximation 

which is accurate to within 0.6 % for < x < 3.273687 {ip{x) diverges at a; sa 3.27368734). This leads to an effective 
charge given by 

This is always less than the actual charge Z which one may interpret as a certain fraction of counterions being 
associated near the surface if A > 0. The effective charge Zeff tends to Z as A ^ (for a general analysis of this 
limit, see the Appendix). 

We now wish to analyze the thermodynamic properties of the crystal in the presence of simple salt which we do 
within a Donnan approximation. At this stage it is well to recall the incorrectness of applying Donnan arguments 
to a fluid of charged colloidal particles. The probability of the double layers of two particles interpenetrating is very 
small owing to Boltzmann weighting. Hence, only the Debye-Hiickel tails in their interaction are important which 
represent effectively the potential of mean force between the particles. In the case of excess salt, we then use the 
McMillan-Mayer theory to calculate the statistical mechanical properties of the fluid as has been done in section II (see 
Eq. (0); this line of argumentation goes back to Stigter [23). The situation is decidedly different when the particles 
are positionally ordered as in a crystal. The double layers are forced to overlap in that case. A usual (Donnan) 
approximation is then to suppose those points at zero electric field are in equilibrium with the reservoir _27J . For the 
cavities in the crystal, this yields 

W + Ic)^!^ (29) 

in view of the equality of the chemical potentials of the small ions in the respective phases. The osmotic pressure is 
given by the additivity rule as argued by Oosawa for polyions within conventional cell models |27l | 

n = {-p + 2Ic - 2T)kB T 

= pkBT[VT+^-w] ^'^^> 

2/ 

w = ^. 
P 

The ions have been considered as ideal and the electrostatic stress is zero in Eq. H20|) . The chemical potential of the 
charged cavity, accurate to the same level of approximation, is readily computed from Eq. H3U|I (this is analogous to 
similar calculations for cell models of long charged rods |30|) 



ZeffksT 
fj. = ^ref H ^ hi 



Vl+w2 + l 



vTTw2 - 1 



(31) 



where firef is a reference chemical potential independent of the concentration of salt, and not identical with /iq of 
section II. Because the number of particles in the crystal is equal to the number of cavities, Eq. (|31|l also represents 
the chemical potential of a charged sphere carrying Z charges but with a different pref ■ 

A. Comparison with experiment 

1. Lysozyme crystal 

The volume per lysozyme molecule in the tetragonal crystal is 29.7 nm'^ (see section III.B). The radius of the 
effective sphere is 1.61 nm so the volume of a cavity is 12.2 nm^^ and b — 1.43 nm. In Tables VIII-X, we show values of 
Z as, a. function of the ionic strength I at three values of the pH. From these we calculate the dimensionless quantities 
A, h/\ and Z^ff via the Poisson-Boltzmann equation. Then the pressure and the chemical potential are evaluated 
using Eqs. (|3()(l and H31|) . (See Tables VIII-X). The curves in Fig. 7 represent the chemical potential computed in 
this manner together with the predictions from the theory of the liquid state as displayed in Fig. 2. The former have 
been shifted by an amount which is unknown in the present theory. 



2. STA crystals 

In order to compute the chemical potential we first need to discuss the crystal habits of STA. It is known that 
H4STA is fully dissociated for a pH larger than 5 [31|. Zukoski et al. [3,|2^|32| fail to mention the pH at which their 
measurements were performed, though they did deduce that all forms of STA are dissociated in their experiments 
judging from the conductivities of their solutions. 

a. H4STA.3IH2O 

This crystallizes at room temperature 133' in the tetragonal form (long axis — 1.856 nm, short axes = 1.301 nm 
|33| : there are 2 STA molecules per unit cell of 3.142 nm^). We have earlier set the radius of an STA ion equal to 
0.55 nm (see section III.C) so the volume of II2O per STA molecule is 0.874 nm'^ or 5 = 0.593 nm. In Refs. p and 
[3^ the water content of this crystal is given in terms of the molecular formula II4STA.3III2O. 

b. LHSTA.24H2 0/LiaSTA . 26H2 O 

Kraus describes two forms of Li4STA with 24 II2O and 26 II2O molecules attached respectively |33- Both crystals 
are rhombohedral (short axes in both cases — 1.559 nm, long axis — 3.898 nm in the former, long axis = 4.118 nm 
in the latter; the angle between the short axes — 120° [23 )■ Kraus also mentions that one Li ion should probably be 
replaced by one H ion. There are actually three different numbers quoted for the water content of Li4STA.nIl2 in 
Refs. [a,|2^|33|: n = 21,24 and 26! We have opted for n ~ 25, namely the average number for the crystal habits 
generally accepted. As there are 6 STA molecules per unit cell, the volume of crystal per STA molecule is 1.406 nm^ 
and the radius of our effective cavity is 6 = 0.553 nm. 

c. NaiSTA.lSHiO 

This crystallizes in the triclinic form within a narrow range around 308 K '33] . The absolute dimensions of the unit 
cell do not seem to be known. We thus estimate the amount of II2O per STA molecule via the molecular formulas. 
The water content in Na4STA.nH2 is stated to be n = 18 in Ref. ;32i and n — 14 in Ref. |^. The latter value seems 
too low and is possibly a misprint since n should be equal to 20 according to the usual citation '36l . Accordingly, we 
adopt ri = 18 here to be used in the solubility studies 25] . A molecule of H2O has a volume of 0.0285 nm"^ which is 
based on the amount of H2O in the unit cells of H4STA.3IH2O and Li4STA.26H20. Therefore, Na4STA.18H20 has 
0.513 nm'^ II2O per STA molecule so we have b = 0.496 nm. 

Overall, it is not clear how much II2O is exactly present in the STA crystals. Fortunately, the chemical potential 
(Eq. (|31|) ') depends only logarithmically on this quantity so the data compiled in Tables XI-XIII are not so sensitive 
to this type of uncertainty. The predicted chemical potentials are depicted as curves in Fig. 8 together with the 
computations from our theory of the liquid state (Fig. 5). 

V. DISCUSSION 

Except for a slight downward adjustment of the effective charge of lysozyme in the fluid phase, there are essentially 
no adjustable parameters in our analysis. The adhesion parameters are completely constrained by the 2nd virial curves 
(Figs. 1 and 4). We predict that the chemical potentials in the fluid and solid phases should coincide apart from 
an unimportant shift in the vertical offset because the reference potential is not known exactly for the crystal. This 
appears to be almost the case for lysozyme (see Fig. 7) but there is an appreciable disparity between the respective 
curves in the case of the silicotungstates (see Fig. 8). Nevertheless, we note that the shapes of the curves are the 
same which implies that the logarithmic form in Eq. (|31(l appears to be confirmed i.e. the Donnan effect seems to 
apply to crystals of charged nanoparticles. This is borne out by adjusting Z^ff upward somewhat for both types of 
crystals. We then actually attain coincident curves (see Figs. 7 and 8). A further implication is that the precise crystal 
structure is unimportant with regard to the ionic-strength dependence of /i. Eq. (|31|l results from approximating the 
cavities within the crystals by spheres; the detailed electrostatics is independent of the salt concentration. 

The coexistence equation for the osmotic pressure yields little information (see Eq. H20|l ') because it is unclear 
how to relate the adhesive forces between spheres to the elastic properties of the crystal. To compute the latter we 



need insight in the forces between the particles at the Angstrom level which we do not have at present. Adhesive 
interactions appear to play a minor role in the STA crystals for the fluid and crystal pressures are quite close (compare 
Tables V-VII with XI- XIII). By contrast, in lysozyme crystals the osmotic pressure due to electrostatic forces is largely 
balanced by sticky interactions between touching protein molecules. 

It is wise to emphasize the shortcomings in the approximations introduced in the electrostatic interactions. Discrete 
charge effects have been disregarded entirely. At the same level of approximation we have not addressed the electro- 
statics of the minutely thin boundary layers between almost touching spheres within the crystal phase. There are 
cavities at nanometer scales and these are assumed to give rise to the ionic-strength dependence of the free energy of 
the crystal. The Donnan approximation used suffers from the same drawback as always: the effective charge .^e// is 
posited to be independent of the electrolyte in the crystal and thus the reservoir (the fluid phase in our case). It would 
be interesting to study the fluid-crystal coexistence of globular particles of low charge density. The counterions in the 
crystal would then be essentially free (see the Appendix) and there would be less uncertainty about the magnitude of 
the electrostatic interactions. 

There is another potential problem in the fluid phases of silicotungstates. At 1 M electrolyte, the solubilities of STA 
are remarkably high (see Tables V-VII) . It would appear that the counterions arising from STA should contribute to 
the screening on a par with the salt ions. This is not borne out by the present analysis, however, since there is no 
levelling off of the chemical potentials in the crystal phases in Figs. 7 and 8. Nevertheless, a liquid state theory of 
concentrated charged nanoparticles needs to be developed in which the counterions are duly accounted for. We note 
that the interaction between the particles is not pairwise additive in that case. 

In summary, we have provided a semi-quantitative explanation for the ionic-strength dependence of the fluid-crystal 
coexistence of suspensions of charged nanoparticles. We believe this explanation is especially forceful because we have 
considered two rather disparate types of globular particles in detail. In particular, the solubility curves of lysozyme 
and silicotungstate differ markedly, yet the curves for the chemical potentials turn out to have the same form. 

VI. APPENDIX: POISSON-BOLTZMANN EQUATION IN A CRYSTAL OR POROUS MEDIUM 

It is possible to present a general analysis of the Poisson-Boltzmann equation for the electrostatic potential ^(r) 
at position r within the aqueous interstitial space inside a crystal (which may be considered to be a porous medium) , 
under appropriate conditions. The particles in the crystal are positively charged and simple salt is absent at first. 
The potential is again related to the countcrion density p[f) via the Poisson Eq. \2.'i\ . Now it is possible to discern 
some point P in the void between several particles where the potential is a local minimum and where the density is 
Pp(0) (see Fig. 9). Point P is chosen as the origin. 

If the potential is scaled analogously as in section IV, we have p{f) = ppexp-0(r) (see Eq. (|22|) '). Thus, the 
Poisson-Boltzmann equation may be written as 

At/- = Ap^e'/' (32) 

where the screening length Xp is given by \p — A'KQ'pp. 

In general, it is difficult to address Eq. (|5^ because Xp is unknown. But it is possible to progress if we suppose 
Ml ^ ^p where \d\ is the largest vector distance between P and a point on the surface of the surrounding spheres (i.e. 
those belonging to a cluster enclosing the void centered on P). A solution of Eq. H32() must have the form ?/'(r/Ap) 
and may be written as a Taylor expansion to second order 



■0(r) — —rr 



1 _ a-^v 



2 drdf 

if \d\ < Xp. Next, we have a boundary condition on the electric field at f— d 



(33) 



^ dip 
or 



fid 



r=d ^^^^ 



= A-KkiCJbQ. (34) 

f=0 



Here, Ub is the uniform density of charge on a sphere and ki is a numerical coefficient of order unity. The effect 
of an internal permittivity is disregarded. The left-hand side of Eq. (|34|l scales as Ap implying that pp must be 
proportional to ajy. In view of electroneutrality we also require the average of p{r) to be proportional to Cf,. Hence, 
the potential ipir) must be very small, which is consistent with the initial Ansatz Eq. (|33|l . We conclude that for 
small enough cavitylike voids, the density of counterions is approximately constant so that the effective charge density 
is virtually equal to the actual charge density. In that case, when the crystal is immersed in a reservoir containing 
monovalent electrolyte, Eqs. H3U|) and H31|) are valid with p simply given by the concentration of counterions in the 
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interstitial space between the spheres; ^e// = ^ in Eq. (|31|l . Because \d\ = 0{a), we ultimately require ZQ/a <C 1 as 
a necessary and sufficient condition for this to hold true. In the spherical cavity approximation introduced in Section 
IV, we have b = 0{a) so A <C 1 is effectively the same requirement (which led to Z^jf = Z). 
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Tables 



/ (%w/v) 


/(M) 


UJ 


Z 


^e// 


Z 


i 


r 


S (g/1) 


V 


knT 


n-uo 

fenT 


2.0 


0.40 


3.35 


11.1 


10.93 


9.93 


1.146 


0.138 


48.7 


0.0358 


-3.55 


0.0319 


3.0 


0.57 


4.00 


11.2 


11.07 


10.07 


0.893 


0.101 


14.0 


0.0103 


-4.69 


0.0097 


4.0 


0.74 


4.56 


11.4 


11.29 


10.29 


0.755 


0.087 


4.30 


0.00316 


-5.80 


0.0031 


5.0 


0.92 


5.05 


11.6 


11.51 


10.51 


0.663 


0.079 


3.11 


0.00229 


-6.12 


0.0022 


7.0 


1.26 


5.92 


11.7 


11.63 


10.63 


0.519 


0.070 


1.36 


0.00100 


-6.93 


0.0010 



TABLE I: The charge Z of hen-egg-white lysozyme (from Rcf. 20]), the effective charge Zeff (from Eq. OJ, the 
lowered effective charge Z = Z^ff — 1, the dimensionless interaction parameters u, ^ and t, the solubility of lysozyme 
S, the volume fraction rj, the dimensionless chemical potential {^ — fj,o)/kBT and the dimensionless pressure Hvo/ksT 
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as a function of the ionic strength / in the fluid phase. The pH equals 4.0 and ^ has been calculated using the lowered 
effective charge Z. 



I (%w/v) 


/(M) 


UJ 


Z 


^e// 


Z 


e 


r 


S (g/1) 


V 




knT 


2.0 


0.40 


3.35 


10.2 


10.07 


9.07 


0.956 


0.114 


30.1 


0.0221 


-4.02 


0.0199 


3.0 


0.57 


4.00 


10.3 


10.20 


9.20 


0.745 


0.089 


10.3 


0.00759 


-4.99 


0.0072 


4.0 


0.74 


4.56 


10.3 


10.22 


9.22 


0.606 


0.078 


5.22 


0.00385 


-5.63 


0.0037 


5.0 


0.92 


5.05 


10.4 


10.33 


9.33 


0.523 


0.072 


3.43 


0.00253 


-6.03 


0.0025 


7.0 


1.26 


5.92 


10.4 


10.35 


9.35 


0.401 


0.065 


1.87 


0.00137 


-6.62 


0.0014 



TABLE II: Same as Table I, but now with a pH equal to 4.5. 



/ (%w/v) 


/(M) 


u; 


Z 


Zeff 


Z 


e 


T 


S (g/1) 


V 




knT 


2.0 


0.40 


3.35 


9.1 


9.00 


8.00 


0.745 


0.094 


17.6 


0.0130 


-4.52 


0.0119 


3.0 


0.57 


4.00 


9.1 


9.03 


8.03 


0.568 


0.077 


7.38 


0.00543 


-5.31 


0.0052 


4.0 


0.74 


4.56 


9.2 


9.14 


8.14 


0.473 


0.070 


4.72 


0.00347 


-5.73 


0.0034 


5.0 


0.92 


5.05 


9.2 


9.15 


8.15 


0.399 


0.066 


3.63 


0.00267 


-5.98 


0.0026 


7.0 


1.26 


5.92 


9.1 


9.07 


8.07 


0.299 


0.060 


2.46 


0.00181 


-6.36 


0.0018 



TABLE III: Same as Table I, but now with a pH equal to 5.4. 



/(M) 


LO 


Z 


^e// 


e 


0.3 


0.99 


4.00 


3.42 


1.905 


1.0 


1.80 


4.00 


3.58 


1.054 


3.0 


3.13 


4.00 


3.76 


0.536 


4.0 


3.61 


4.00 


3.80 


0.439 


5.0 


4.03 


4.00 


3.83 


0.373 



TABLE IV: Values of the bare charge Z of STA, the effective charge ^e// (from Eq. ^) and the dimensionless 
interaction parameters to = 1.80\/T and ^ — 0.645(Zey//(l + cj))^ as a function of the ionic strength /. These entries 
apply to the data plotted in Fig. 4. 



/(M) 


U) 


Z 


Zeff 


e 


T 


S (g/ml) 


V 


knT 


knT 


1.0 


1.80 


4.0 


3.58 


1.054 


0.786 


1.94 


0.284 


1.30 


0.706 


2.0 


2.55 


4.0 


3.70 


0.700 


0.358 


1.67 


0.243 


-0.48 


0.375 


3.0 


3.13 


4.0 


3.76 


0.536 


0.255 


1.36 


0.198 


-1.40 


0.226 


4.0 


3.61 


4.0 


3.80 


0.439 


0.215 


1.11 


0.162 


-1.91 


0.158 


5.0 


4.03 


4.0 


3.83 


0.373 


0.193 


0.57 


0.0830 


-2.65 


0.077 



TABLE V: The charge Z of STA, the effective charge Z^ff (from Eq. (j2Jl), the dimensionless interaction parameters 
w, f and T, the solubility S of II4STA, the volume fraction 77, the dimensionless chemical potential (/i — ^oj/ksT and 
the dimensionless pressure nwo/Zc^T as a function of the ionic strength I in the fluid phase. Here the counterion is 
H+ and the added salt is HCl. ^ has been calculated using the effective charge Z^ff- 



/(M) 


LO 


Z 


^e// 


i 


r 


S (g/ml) 


V 


knT 


ni;o 
knT 


1.0 


1.80 


4.00 


3.58 


1.054 


0.668 


2.15 


0.312 


1.55 


0.811 


2.0 


2.55 


4.00 


3.70 


0.700 


0.352 


1.85 


0.267 


-0.26 


0.433 


3.0 


3.13 


4.00 


3.76 


0.536 


0.255 


1.36 


0.197 


-1.41 


0.224 


4.0 


3.61 


4.00 


3.80 


0.439 


0.215 


0.81 


0.118 


-2.23 


0.114 


5.0 


4.03 


4.00 


3.83 


0.373 


0.193 


0.32 


0.0469 


-3.16 


0.045 
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TABLE VI: Same as Table V, but now with Li"*" as counterion and LiCl as added salt. 



I (M) 


UJ 


Z 


^e// 


C 


T 


S (g/ml) 


'7 


M-Mo 
k„T 


nuo 

k„T 


1.0 


1.80 


4.0 


3.58 


1.054 


0.906 


1.85 


0.262 


1.07 


0.625 


2.0 


2.55 


4.0 


3.70 


0.700 


0.359 


1.67 


0.236 


-0.54 


0.358 


3.0 


3.13 


4.0 


3.76 


0.536 


0.255 


1.36 


0.192 


-1.44 


0.218 


4.0 


3.61 


4.0 


3.80 


0.439 


0.215 


0.74 


0.104 


-2.35 


0.100 


5.0 


4.03 


4.0 


3.83 


0.373 


0.193 


0.46 


0.0657 


-2.86 


0.062 



TABLE VIL Same as Table V, but now with Na+ as counterion and NaCl as added salt. 



/ (%w/v) 


/(M) 


Z 


A 


h/X 


Zeff 


w 


k„T 


knT 


2.0 


0.40 


11.1 


5.51 


2.52 


4.27 


1.39 


2.87 


1.97 


3.0 


0.57 


11.2 


5.56 


2.53 


4.29 


1.97 


2.10 


1.47 


4.0 


0.74 


11.4 


5.66 


2.54 


4.32 


2.54 


1.66 


1.17 


5.0 


0.92 


11.6 


5.76 


2.55 


4.35 


3.10 


1.38 


0.98 


7.0 


1.26 


11.7 


5.81 


2.55 


4.36 


4.24 


1.02 


0.72 



TABLE VIIL The ionic strength /, the actual number of charges Z, the effective number ^e//, the chemical 
potential /i and the osmotic pressure H for a lysozyme crystal at pH 4.0 (the reference chemical potential has been 
set equal to zero). 



/ (%w/v) 



2.0 



3.0 



4.0 



5.0 



7.0 



/(M) 



0.40 



0.57 



0.74 



0.92 



10.2 



10.3 



10.3 



10.4 



1.26 10.4 5.17 



A 



5.07 



5.12 



5.12 



5.17 



6/A 



2.48 



2.49 



2.49 



2.49 



2.49 



eff 



4.13 



4.15 



4.15 



4.16 



4.16 



1.43 



2.04 



2.64 



3.24 



4.45 



fesT 



2.69 



1.96 



1.53 



1.27 



0.93 



TABLE IX: Same as Table VIII, but now for pH 4.5. 



TABLE X: Same as Table VIII, but now for pH 5.4. 



keT 



1.86 



1.38 



1.08 



0.90 



0.66 



/ (%w/v) 


/(M) 


Z 


A 


b/\ 


Zeff 


w 


fesT 


knT 


2.0 


0.40 


9.1 


4.52 


2.42 


3.93 


1.51 


2.45 


1.70 


3.0 


0.57 


9.1 


4.52 


2.42 


3.93 


2.15 


1.77 


1.25 


4.0 


0.74 


9.2 


4.57 


2.43 


3.95 


2.77 


1.40 


0.99 


5.0 


0.92 


9.2 


4.57 


2.43 


3.95 


3.41 


1.14 


0.81 


7.0 


1.26 


9.1 


4.52 


2.42 


3.93 


4.71 


0.83 


0.59 



/ (%w/v) 


/(M) 


Z 


A 


b/X 


Zeff 


w 


k„T 


knT 


1.0 


0.60 


4.0 


4.789 


2.452 


1.67 


0.629 


2.08 


0.738 


2.0 


1.20 


4.0 


4.789 


2.452 


1.67 


1.257 


1.22 


0.466 


3.0 


1.81 


4.0 


4.789 


2.452 


1.67 


1.886 


0.85 


0.332 


4.0 


2.41 


4.0 


4.789 


2.452 


1.67 


2.514 


0.65 


0.256 


5.0 


3.01 


4.0 


4.789 


2.452 


1.67 


3.143 


0.52 


0.207 



TABLE XI: Same as Table VIII but now for H4STA.3IH2O. 
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/ (%w/v) 


/(M) 


Z 


A 


6/A 


^e// 


w 




Wvn 
k„T 


1.0 


0.60 


4.0 


5.134 


2.488 


1.61 


0.531 


2.23 


0.950 


2.0 


1.20 


4.0 


5.134 


2.488 


1.61 


1.062 


1.35 


0.626 


3.0 


1.81 


4.0 


5.134 


2.488 


1.61 


1.594 


0.95 


0.455 


4.0 


2.41 


4.0 


5.134 


2.488 


1.61 


2.125 


0.73 


0.353 


5.0 


3.01 


4.0 


5.134 


2.488 


1.61 


2.656 


0.59 


0.288 



TABLE XII: Same as Table VIII but now for Li4STA.25H20. 



I (%w/v) 


/(M) 


Z 


A 


b/X 


Zeff 


w 


k„T 


ni>o 

knT 


1.0 


0.60 


4.0 


5.721 


2.542 


1.51 


0.410 


2.45 


1.374 


2.0 


1.20 


4.0 


5.721 


2.542 


1.51 


0.819 


1.55 


0.970 


3.0 


1.81 


4.0 


5.721 


2.542 


1.51 


1.229 


1.12 


0.728 


4.0 


2.41 


4.0 


5.721 


2.542 


1.51 


1.639 


0.87 


0.575 


5.0 


3.01 


4.0 


5.721 


2.542 


1.51 


2.049 


0.71 


0.473 



TABLE XIII: Same as Table VIII but now for Na4STA.18H20. 



Figure Captions 



FIG. 1. The second virial coefficient of lysozyme as a function of the ionic strength. The second virial coefficient is 
scaled by the hard sphere value 82^ . The data are taken from a variety of experiments; see Ref. T\ for more details. 
The added salt is NaCl. The solid hne is a fit to the data with Ua = 2.89 and 5 = 0.182. 

FIG. 2. The dimensionless chemical potential of lysozyme in the fluid phase as a function of the ionic strength at 
pH 4.0 (diamonds), pH 4.5 (squares) and pH 5.4 (triangles). See also Tables I, II and III. 

FIG. 3. The dimensionless pressure of lysozyme in the fluid phase as a function of the ionic strength at pH 4.0 
(diamonds), pH 4.5 (squares) and pH 5.4 (triangles). See also Tables I, II and III. 

FIG. 4. The second virial coefficient of STA as a function of the ionic strength. The second virial coefficient is 
scaled by the hard sphere value B2^. The experimental data are taken from Zukoski et al. (Ref. [23 )• The added 
salt is LiCl (diamonds), HCl (squares) and NaCl (triangles) respectively and in all cases the counterion of STA is the 
same as that of the salt. The solid line is a fit to the experimental data with Ua = 3.30 and 5 = 0.05. 

FIG. 5. The dimensionless chemical potential of STA in the fluid phase as a function of the ionic strength. The 
added salt is LiCl (diamonds) , HCl (squares) and NaCl (triangles) respectively and in all cases the counterion of STA 
is the same as that of the salt. See also Tables V, VI and VII. 

FIG. 6. The dimensionless pressure of STA in the fluid phase as a function of the ionic strength. The added salt is 
LiCl (diamonds), HCl (squares) and NaCl (triangles) respectively and in all cases the counterion of STA is the same 
as that of the salt. See also Tables V, VI and VII. 

FIG. 7. The chemical potential of lysozyme in the fluid phase as a function of the ionic strength at pH 4.0 
(diamonds), pH 4.5 (squares) and pH 5.4 (triangles) (see Fig. 2). The solid lines denote predictions from the theory 
of the crystalline state (Eq. IJ3U), with the effective charge from Table VIII-X. The shift in chemical potential in 
units of UbT has been chosen to be 7.4 (light grey line, pH 4.0), 7.3 (black line, pH 4.5) and 7.15 (dark grey line, pH 
5.4) respectively. Dashed line denotes the theory of the crystal for Z^ff — 5.0 (shift — 7.9) and the dash-dotted line 
for Zeff = 5.9 (shift = 8.7). 

FIG. 8. Chemical potential of STA in the fluid phase as a function of the ionic strength (see Fig. 5). Salt added: 
LiCl (diamonds); HCl (squares); NaCl (triangles) (counterion of STA is the same as that of the salt). Solid lines 
denote predictions from the theory of the crystalline state (Eq. IJ3U), with the effective charge from Table XI-XIII. 
The shift in chemical potential in units of fc^T is 2.3 (black hne, H4STA.3IH2O), 2.4 (dark grey hne, Li4STA.25H20) 
and 2.55 (light grey line, Na4STA.18H20) respectively. Dashed line denotes predictions from the theory of the crystal 
for Zefs = 2.8 (shift = 5.5). 

FIG. 9. Point P in a void of the crystal. 
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